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A SQUARE ROOT EXTRACTOR 

Technical Field of the Invention 

The present invention relates in general to an 
arrangement for producing an output signal which is the 
5 square root of an input signal . 



Background of the Invention 

A wide variety of applications require the use 
of a square root extractor. For example, at least some 
electric-force-balance instruments are non-linear devices 

10 used to measure various physical phenomena. Such 

instruments typically develop an internal restoring force 
which is proportional to the square of an applied analog 
control signal. The analog control signal is usually 
provided by a digital controller, which supplies a 

15 digital command signal, and a digital-to-analog converter 
(DAC) , which converts the digital command signal to an 
analog command signal . This analog command signal is 
supplied to a driver which provides the analog control 
signal to the electric-force-balance instrument. 

2 0 As suggested by the above, the response of the 

electric-force-balance instrument varies as the square of 
the amplitude of the control signal. Accordingly, the 
response of the electric- force-balance instrument is non- 
linear with respect to its input control signal. Because 

2 5 of this non-linearity, linear control algorithms cannot 
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be used directly to control the electric-force-balance 
instrument. Instead, such linear control algorithms must 
incorporate additional processing which linearizes the 
response of the instrument, thereby adding extra cost and 
5 complexity to the electric-force-balance instrument. 

The present invention, therefore, is directed 
to an arrangement for providing an output signal that has 
an amplitude which is proportional to the square-root of 
the amplitude of an input signal. This arrangement is 
10 advantageous because it permits the direct use of linear 
control algorithms to control instruments having non- 
linear responses. 



Summary of the Invention 

In accordance with one aspect of the present 

15 invention, an apparatus comprises a sign extractor, a 

delay, a square root extractor, and a sign restorer. The 
sign extractor has an input and first and second outputs. 
The input of the sign extractor receives an input signal, 
the first output provides a sign of the input signal, and 

2 0 the second output provides a magnitude of the input 

signal. The delay is coupled to the first output, and 
the delay imposes a delay on the sign. The square root 
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extractor is coupled to the second output . The square 
root extractor has an output that provides an output 
signal, and the output signal is an approximation to a 
square root of the magnitude of the input signal. The 
5 sign restorer is coupled to the output of the square root 
extractor and to the delay. The sign restorer applies 
the sign from the delay to the output signal from the 
square root extractor. 

In accordance with another aspect of the 

10 present invention, a square root extractor consists of 
multiplying, summing, scaling, and delaying functions. 

In accordance with yet another aspect of the 
present invention, a square root extractor comprises 
first, second, and third multipliers, first and second 

15 summers, first and second delays, and a scaler. The 

first multiplier has a first input coupled to receive a 
signal whose square root is to be extracted and a second 
input coupled to an output of the third multiplier. The 
first summer has a first input coupled to an output of 

20 the first multiplier and a second input coupled to an 
output of the first delay. The scaler has an input 
coupled to an output of the first summer and an output 
coupled to an input of the first delay. The output of 
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the first delay provides an output of the square root - 
extractor. The second multiplier has a first input 
coupled to the output of the first delay and a second 
input coupled to the output of the second delay. The 
5 second summer has a first input coupled to an output of 
the second multiplier and a second input coupled to a 
constant. The third multiplier has a first input coupled 
to an output of the second summer and a second input 
coupled to the output of the second delay. The second 
10 delay has an input coupled to the output of the third 
multiplier. 

In accordance with still another aspect of the 
present invention, a method comprises the following: 
multiplying first and second signals to produce a third 

15 signal, wherein the first signal is a signal whose square 
root is to be extracted; summing the third signal and a 
fourth signal to produce a fifth signal; scaling the 
fifth signal to produce a sixth signal; delaying the 
sixth signal to produce the fourth signal; multiplying 

20 the fourth signal and a seventh signal to produce an 
eighth signal; subtracting the eighth signal from a 
constant to produce a ninth signal; multiplying the 
ninth signal and the seventh signal to produce the second 
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signal; and, delaying the second signal to produce the 
seventh signal, wherein both the fourth signal and the 
sixth signal are approximations to the square root of the 
first signal, and wherein both the second signal and the 
5 seventh signal are approximations to the reciprocal of 
the square root of the first signal. 



Brief Description of the Drawings 

These and other features and advantages will 
become more apparent from a detailed consideration of the 
10 invention when taken in conjunction with the drawings in 
which : 

Figure 1 shows an exemplary operating 
environment for the square root extractor of the present 
invention; and, 
15 Figure 2 shows a square root extractor 

according to one embodiment of the present invention. 



Detailed Description 

As shown in Figure 1, a linear digital 
controller 10 generates a signal on an output 12. As is 
20 typical, the signal on the output 12 from the linear 

digital controller 10 may be used to drive an instrument 
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whose response is proportional to this signal. Thus, 
after digital-to-analog conversion by a digital-to-analog 
converter 14, this signal can be used directly in many 
typical linear-control applications. However, if the 
5 signal on the output 12 is to be used to drive an 

instrument 16 whose internal physics provide a responsive 
force within the instrument that is proportional to the 
square of the magnitude of the signal on the output 12, a 
square root extractor 18 is interposed between the linear 

10 digital controller 10 and the digital-to-analog converter 
14, as shown in Figure 1. The square root extractor 18 
produces a signal on its output 2 0 that is an 
approximation of the square root of the magnitude of the 
signal on the output 12 and carrying the sign of the 

15 signal on the output 12 . 

The square root extractor 18 may be implemented 
in accordance with the following mathematical analysis. 
The magnitude of the signal on the output 12 of the 
linear digital controller 10 may be defined as jx n | and an 

2 0 approximation of its square root may be defined as u n . 

The choice of an approximation error function determines 
the mechanization complexity and the speed of convergence 
provided by the square root extractor 18. However, to 
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simplify mechanization complexity and increase the speed 
of convergence, the selected approximation error function 
may be defined in accordance with the following equation: 



5 The derivative of f (u n ) is then determined in accordance 
with the following equation: 



The zero value for the approximation error function f (u n ) 
may be computed iteratively using the well known Newton- 
10 Raphson update formula from elementary calculus. The 
Newton-Raphson update formula for u is given by the 
following equation: 



f(u n )=ul- \x n \ 



(1) 



/'("„) = 2u n 



(2) 



- u n - 



f{u n ) 



(3) 
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Substituting equations (1) and (2) into equation (3) 
produces the following equation: 




In order to avoid a division step, the method disclosed 
5 in K. Martin, "Power-Normalized Update Algorithm for 
Adaptive Filters Without Division," IEEE ASSP Trans . , 
vol. 37, no. 11, Nov. 1989; pp. 1782-1786 may be used. 
According to this method, l/u n can be approximated as 
v n+1 . Substituting v n+1 for l/u n in equation (4) produces 
10 the following equation: 



: 2\ U " 



The approximation error function for v n may be selected 
in accordance with the following equation: 

1 
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The derivative of equation (6) is given by the following 
equation : 

g'(vj=vf (7) 



The Newton- Raphson update formula for v n is given by the 
following equation: 



g(vj 

g(y n ) 



Substituting equations (6) and (7) into equation (8) 
produces the following equation: 

V «+l = V n{ 2 ~ U n V n) (9) 



10 It should be noted that equation (9) is the reciprocator 
according to the K. Martin paper disclosed above. 

Equations (5) and (9) may be implemented as a 
square root extractor requiring only the simple 
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operations of adding, subtracting, multiplying, and 
delaying in order to produce an output signal that is the 
square root of an input signal. Therefore, a square root 
extractor according to equations (5) and (9) permits the 
5 direct use of linear control algorithms to measure 

observed phenomenon in a simple, straight forward manner. 

Because the operation required by equations (5) 
and (9) is iterative, there are practical bandwidth 
restrictions in using a square root extractor according 

10 to these equations. However, from arbitrary initial 

conditions (such as u 0 = 0 and v 0 = 0.001), convergence 
between u n , the square root approximation of |x n | , and the 
actual square root of \x n | within parts per billion is 
achieved in less than a dozen sample periods. When 

15 following a well-behaved signal such as a sinusoid, the 
tracking error is very small. 

Figure 2 shows an implementation of the square 
root extractor 18 in accordance with equations (5) and 
(9) . The signal (XjJ on the output 12 of the linear 

20 digital controller 10 is coupled to a sign extractor 40 

having first and second outputs 42 and 44. The signal on 
the first output 42 of the sign extractor 40 is the sign 
of the signal (x n ) on the output 12 . The signal on the 
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second output 44 of the sign extractor 4 0 is the 
magnitude of the signal (xj on the output 12. That is, 
the signal on the second output 44 of the sign extractor 
4 0 is the absolute value of the signal (x n ) on the output 
5 12 . In this regard, the sign extractor 40 may be 

arranged to complement the signal (x n ) and to provide 
either the signal (x n ) or the complement of the signal 
(Xn) on the second output 44 depending on which of the 
signals (x n ) has the positive sign bit. 

10 The second output 44 is coupled to a first 

input of a first multiplier 46. The signal from the 
output of the first multiplier 46 is summed by a first 
summer 48 with a signal produced by a first one-sample- 
period-delay element 50. The output from the first 

15 summer 48 is scaled by y 2 by a scaler 52 . The output of 
the scaler 52 is coupled to an input of the first one- 
sample-period-delay element 50, and the output of the 
first one- sample-period-delay element 50 is an 
approximation of the square root of the magnitude of the 

20 signal (x n ) on the output 12 . The output of the first 

one-sample-period-delay element 50 is provided to a first 
input of a sign restorer 54. The sign on the output 42 
from the sign extractor 40 is delayed by a second one- 
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sample-period-delay element 56 and is coupled to a second 
input of the sign restorer 54. The sign restorer 54 
merely applies the sign from the second one-sample- 
period-delay element 56 to the output signal at the 
output of the first one- sample-period-delay element 50. 
Thus, the signal provided by the sign restorer 54 on the 
output 20 is an approximation of the square root of the 
magnitude of the amplitude of the signal (x n ) on the 
output 12 and has the sign of the signal (x n ) on the 
output 12 . 

The output from the first one -sample -period- 
delay element 50 is further coupled to a first input of a 
second multiplier 58. The second multiplier 58 produces 
an output signal which is coupled to a negative input of 
a second summer 60. A constant k = 2 is provided to a 
positive input of the second summer 60. The second 
summer 60, accordingly, subtracts the output of the 
multiplier 58 from the constant k = 2. The second summer 
60 produces an output signal which is coupled to a first 
input of a third multiplier 62. The third multiplier 62 
produces an output signal which is coupled to a second 
input of the first multiplier 46. The output signal from 
the third multiplier 62 is also delayed by a third one- 
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sample-period-delay element 64. The third one-sample- 
period-delay element 64 provides a signal on an output 66 
which is an approximation of the reciprocal of the square 
root of the signal at the second output 44 from the sign 
5 extractor 40. The output 66 is coupled to second inputs 
of the second multiplier 58 and the third multiplier 62. 

Accordingly, the square root extractor 18 
implements equations (5) and (9) to produce an 
approximation of the square root of the magnitude of the 

10 signal on the output 12 from the linear digital 
controller 10. 

Certain modifications of the present invention 
will occur to those practicing in the art of the present 
invention. For example, the square root extractor 18 can 

15 be a digital square root extractor, as shown and 

described above, or the square root extractor 18 can be 
an analog square root extractor. If the square root 
extractor 18 is an analog square root extractor, it may 
be desirable to interpose the square root extractor 18 

20 between the digital-to-analog converter 14 and the 

instrument 16. Alternatively, the square root extractor 
18 having an analog form may be used without the digital- 
to-analog converter 14 in the case where a linear analog 
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controller is used in place of the linear digital 
controller 10. 

Accordingly, the description of the present 
invention is to be construed as illustrative only and is 
for the purpose of teaching those skilled in the art the 
best mode of carrying out the invention. The details may 
be varied substantially without departing from the spirit 
of the invention, and the exclusive use of all 
modifications which are within the scope of the appended 
claims is reserved. 
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